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ABSTRACT 

A fine tuned hierarchy between a strongly coupled high energy compos- 
iteness scale and a much lower chiral symmetry breaking scale is a requisite 
ingredient in many models of dynamical electroweak symmetry breaking. 
Using a nonperturbative continuous Wilson renormalization group equation 
approach, we explore the stability of such a hierarchy against quantum fluc- 
tuations. 



Many of the currently studied models of dynamical electroweak symmetry 
breaking include some strongly interacting sector acting at a high energy 
scale, A > lOTe^, which produces an essentially composite scalar bosonic 
degree of freedom. In addition, this dynamics is also supposed to play a non- 
trivial role in the electroweak symmetry breaking whose characteristic scale 
is much lower; A^? ~ 250GeV^. Thus these models require that a significant 
hierarchy can be established between these scales. Included in this class are 
strong extended technicolor models'^', models of broken technicolort^l and 
models involving heavy quark condensation^^"^!. 
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In general, the hierarchy is achieved by a fine tuning of parameters close 
to the critical value for the chiral symmetry breaking. A prototype of this 
behavior is exhibited by the Nambu Jona-Lasinio (NJL) model 1^1, where a 
fine tuning of the four-fcrmion coupling allows the emergence of a chiral 
symmetry breaking scale far below the compositeness scale. A feature of this 
particular model which must be generic in any model which exhibits such a 
large hiearchy of scales is that the chiral symmetry phase transition be of 
second order. That is, in order for the hierarchy to be maintained and not 
have the electroweak scale driven to be of order A, it is necessary that the 
order parameter characterizing the chiral transition must remain zero as the 
theory is scaled from A into the infrared until one reaches the electroweak 
scale. 

If, on the other hand, quantum fiuctuations turn the transition first order 
at a scale e~*°A >> A^, then the order parameter will jump discontinuously 
to be of this value and it will be impossible to maintain the hierarchy all 
the way down to the electroweak scale. Instead the hierarchy will destabilize 
after to e-foldings. Such a situation is an example of the phenomenon studied 
by Coleman and Weinbergt^l . It is important to recognize that this question 
is distinct from that of the naturalness of the fine tuning of additive quadratic 
divergences. Even allowing for such fine tunings, this destabilization of hi- 
erarchies could prove problematic for many of the models of dynamical elec- 
troweak symmetry breaking. Thus, as emphasized by Chivukula, Golden and 
Simmons'^l, it becomes essential to explore when the hierarchy can be self 
consistently maintained and not destroyed by quantum fiuctuations. 

The minimal NJL model which is used in the minimal top condensate 
model'^1 is known to exhibit a second order chiral transition. However, since 
a top quark mass of 174 GeV, as reported by the CDF coUaborationl^^l, is 
outside the range of the minimal model, it appears that if such models are 
to be phenomenologically viable, they must necessarily require extensions 
beyond the minimal version. Similarly, while the 0(4) linear sigma model 
containing a single scalar quartic self coupling has been shown to exhibit 
a mean field second order chiral transition, the effective dynamics of strong 
extended technicolor models could require an effective Lagrangian containing 
more than one scalar self couphng. Thus we are led to investigate models 
containing multiple scalar quartic self couplings. 

Following previous work I^'^^'i^]^ -^g focus on a model possessing a global 
chiral U{2)lxU (2)ij symmetry which has two independent scalar quartic self 
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couplings. The model degrees of freedom include left and right handed chiral 
fermions ipii and ipm, i = 1,2, transforming as the fundamental (2, 0) and 
(0, 2) representations of the [7(2)^ x C/(2)r groups respectively which further 
carry the fundamental, Nq, representation of an asymptotically free gauged 
symmetry. Throughout our discussion, we shall neglect these gauge interac- 
tions which are considered to be very feeble at these high energy scales. We 
assume that, as a consequence of some unspecified dynamics acting at scale 
A, the U{2)lxU (2)r chiral symmetry is spontaneously broken. This symme- 
try breaking is further assumed to produce a gauge singlet scalar composite 
Hij which has the U{2)l x U{2)r quantum numbers of the fermion bilinear 
ipjRipiL- That is Sjj transforms as the (2,2) under the chiral group. Its vac- 
uum expectation value, can be interpreted as an order parameter for 
the chiral symmetry breaking. Since we are assuming the chiral symmetry 
phase transition is second order, we are led to study a Ginzburg-Landau 
effective model whose Euclidean action at scale A is given by 

5[S,St,^,^;0] = J d'^x{tr{df,E'^d^'^) ■ Dilj + V{x,y,0) 

where Dfj^ip is the fermion gauge covariant derivative and the invariant po- 
tential function takes the form 

V{x, y, 0) = ^m'{0)x + ^XM^' + Y^2{0)y . (2) 

Here x = ir(EtE),y = tr{i:^T,f are the two independent U{2)l x U{2)r in- 
variants. A Coleman- Weinberg instability is signalled by the appearance of a 
non-trivial global minimum of the effective potential, Veff{v), with vanishing 
renormalized mass. If this occurs at the scale v — e~*°A, then the phase 
transition is driven first order by quantum fiuctuations at this scale and one 
can technically achieve a hierarchy of only to e-foldings. 

The conditions for the appearance of such a non-trivial global minimum 
of the effective potential with vanishing renormalized mass can be elegantly 
expressed in terms of the various renormalization group functions by employ- 
ing a formalism due to Yamagishil^^l In general, the (Yamagishi) stability 
function is defined as 
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with V = e The effective potential is minimized provided 

dY 

n^o) = ; ^It=to<0 (4) 

and this minimum is a global one provided V^//(f ) < ^//(O) = 0. Approxi- 
mating the full effective potential by the 1-loop perturbation theory improved 
version which includes only the leading logarithmic radiative corrections of 
the dimension four operators appearing in the potential function of Eq. (2), 
while further holding the Yukawa coupling fixed, these Yamagishi conditions 
for the appearance of a global minimum at to reduce to 

y(to) = [4Ai+4A2 + /3i(A)+/32(A)]|t=to = 

[4 + /^i^ + /^2T-](A + /32)|t=to>0 ; A2(to)>0 (5) 

Ai(io) + A2(io) < 0, 

The first two hnes of Eq. (5) are the fixed Yukawa coupling, 1-loop approx- 
imation version of Eq. (4) and guarentee a minimum at t^. while the last 
inequality insures the minimum is a global one. Note that for other ap- 
proximation schemes which are still restricted to include only the marginal, 
dimension four operators, the Yamagishi coditions of Eq. (5) are modified 
by the replacements — > ^iii) — xtilt) ^^"^ ~^ K{t), where ^{t) is 
the anomalous dimension of E and Pi — Using the 1-loop perturbative 
f3 functions, the Yamagishi stability condition simply reduces to a curve in 
the A2(t)/(y'^(0) — Ai(t)/(y'^(0) plane. The signal for the Coleman- Weinberg 
instability is that the renormalization group flow of the running couplings 
cross this stability curve. Using the 1-loop approximation with fixed Yukawa 
couphng, Chivukula et al.'^' found that, depending on the initial choice of 
couplings, the trajectories either 

(i) run to the infrared quasi fixed point near the origin, 

(ii) cross the stability line signaling a first order transition 

(iii) or simply run away in which case the model is ill defined. 

For example, focusing on the region of coupling space corresponding to 
the initial couplings A2(0) and g^{0) both large and positive and Ai(0) = 
0, the transition goes first order for sufficiently large X2{0)/ g'^{0) ratio 
7). Moreover, this occurs, in general, near to the compositeness scale. For 
example, choosing A2(0) = 10 and g^{0) = 1 (and taking Nc = 3), we display 
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Figure 1: Y{t) as a function of t for the initial parameters Ai(0) = 0, A2(0) = 
10,(7^(0) = 1 computed using the 1-loop renormahzation group improved 
effective potential with fixed Yukawa coupling. 

in Fig. 1 the form of Y{t) as a function of t. This function is seen to cross the 
stabihty curve y = at ~ 1.3. Thus the transition turns first order and 
the hierarchy destabilizes after ~ 1.3 e-foldings (which corresponds to f ~ 
0.27A). However, since the couplings are very large, the 1-loop perturbative 
approximation can certainly be called into question. For instance, naively 
using the 2-loop perturbative renormahzation group functions for these large 
initial couplings, the renormahzation group trajectories do not cross the line 
but simply run away. 

Clearly, some nonperturbative approximation scheme is required to prop- 
erly deal with the system in the vicinity of the strong coupling compositeness 
scale. The purely bosonic U{2)lxU (2)r model (no chiral fermions) has been 
simulated using lattice Monte Carlo techniques^^^l and was seen to undergo 
a Coleman- Weinberg instability. An alternate approach which includes the 
chiral fermions has been advocated by Bardeen et al.'^^^. Modelling the non- 
perturbative physics in the vicinity of the compositeness scale using a large 
Nc approximation, then at scale A the Yukawa coupling is seen to dominate. 
Retaining only it and the fine tuned scalar mass term needed to cancel the 
additive quadratic divergence, the model at scale A reduces to the minimal 
NJL model which is exactly soluble in the large Nc limit producing the run- 
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ning couplings Ai(t) = and A2(t) = 35'^(t) = This model exhibits a 
(trivial) second order chiral transition and as such allows any sized hierarchy 
to be technically achieved. Note that for couphngs consistent with the large 
Nc approximation, the 1-loop perturbative solution is driven to the infrared 
quasi fixed point and thus also does not cross the stability curve. As such it 
too allows for arbitrary large hierarchies to be tuned. Running the couplings 
using the large Nq approximation solution until the couplings have decreased 
sufficiently to be smoothly joined onto a perturbative running which also in- 
cludes the running of the Yukawa coupling, it is found that the stability line 
is eventually crossed but only after a sizeable hierarchy (to ~ 20 — 25) has 
been established. Note that for this case of initially dominate Yukawa cou- 
pling (or for one of comparable size to the initial scalar self coupling), it is 
important to include the effects of its running which is to enhance the flow 
of the system toward the stability line. For smaller initial Yukawa couplings 
the effects of its running will not be nearly as important. Since the large A^^^ 
approximation is nonpcrturbativc, the procedure of Bardeen et al. is a self 
consistent one. On the other hand, it can be reliably employed for only a 
very limited range of the initial parameter space. 

An alternate nonperturbative method is provided by the continuous Wil- 
son renormalization group equation (WRGE)[^^-2o] ^^Ach has been extended 
to include chiral fermions'^^^ . The WRGE nonperturbativcly relates the form 
of the Euclidean action at a scale e~*A to the action at scale A for t > 0. It 
is derived by demanding that the physics, ie. correlation functions, remain 
unchanged as the degrees of freedom carrying momentum between scales A 
and e^*A are integrated out. Thus either action can be used to equivalently 
describe the physics on all scales less than e"*A and both actions lie on the 
same Wilson renormalization group trajectory. The lower scale action is con- 
structed by appropriately changing the coefficients of the operators already 
present at scale A, as well as including new ones, in such a way so as to 
keep the physics unchanged. In general, the resultant action incorporates 
a complete set of local operators. This includes not only the relevant and 
marginal operators, but also irrelevant ones. The coefficient of each operator 
is fixed in terms of the initial action defined at scale A. The full WRGE is 
a very complicated integro functional differential equation and its analysis 
necessarily requires some simplifying approximations. One commonly used 
approximation is to work within a local action approximation^^^^^^l which 
ignores anomalous dimensions and derivative interactions while we further 
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neglect operators higher than bihnear in the fermion fields. In addition, in 
order to obtain a tractable analysis, we are forced to restrict attention to the 
case of a fixed Yukawa coupling. As such, the initial parameter space that we 
are able to investigate using this approach is restricted to be one where the 
initial scalar self coupling(s) dominates the initial Yukawa coupling which, 
however, can still be of order unity. Thus the non-perturbative WRGE we 
solve probes a region of initial parameter space intermediate between the 
lattice Monte Carlo simulations {g{0) = 0) and the large Nc analysis {g{0) 
dominates) and our results can be be viewed as complementary to those 
analyses. Taking into account the various simplifications, the action at scale 
e~*A can be written as 



S[J:, ijj, -0; = / d'^x{tr(d^J:''d^J:) +jjj-fD'ijj + V{x, y, t) 



(6) 



where V{x,y,t) is the U{2)l x U{2)ji invariant potential function at this 
scale. The full WRGE!^'''! then reduces to a partial differential equation for 
this function of the form 



dV 

'at 



dV 



+ 



-,tr en{l + W) 



Nr. 



167r2 ' ' 47r2 
where VF(E, E''') is the matrix of second derivatives 

W 



tr M1 + y/(0)S^S), 



(7) 



d'Eijd'Eke 



aEl,5E^ 



11^ ^kl 

dUijdTi],^ 



(8) 



The tr in terms containing W arise from integrating out the scalar modes, 
while those containing g^{0) are due to integrating out the fermion modes. 
After a considerable amount of algebraic manipulation, the various determi- 
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nants can be explicitly evaluated yielding the WRGEl^^l 



+ + 2x(l + + 2(x2 - 

+ K)(l + K + 6xy^ + 2xKx + SyKy + 4x(3y - x'')Vyy) 
+6(a;2 - |/)K,(3K, + 2V.,^ + 4a;Ky + ^yVyy) 

+ ^x/(0) + ^(x^ - y)^^(O)] , 



47r2 ^ 2 

(9) 

subject to the initial condition of Eq.(2). Here the subscripts denote differ- 
entiation with respect to that variable so that, for example, = etc. For 
t > 0, each action constructed using the V{x^ t) satisfying this equation lies 
on the same Wilson renormalization group trajectory and produces the same 
physics on all scales less than e'~*A. Unfortunately, the solution to this equa- 
tion is currently beyond our numerical abilities. Thus we make the further 
truncation of retaining terms only up to linear in y with coefficients which 
are arbitrary functions of x. Eq. (9) then reduces to two coupled equations 
which are of a similiar (although considerably more complicated) form to 
what we previously solved in obtaining nonperturbative mass bounds'^^'^^^ . 
While the truncations used are drastic and uncontrolled, they still include 
contributions from an infinite number of operators. 

The resulting equations are then numerically solved for t values up to 
some t*, where t* lies in a region where V{x,y,t) is found to be linear in 
t — t* with a slope of the same form as the linearized m. t — t* full 1-loop 
effective potential. The radiative corrections arising from the momentum 
modes less than e^**A can then be satisfactorially incorporated using the 1- 
loop approximation to the effective potential. Note that this 1-loop effective 
potential also includes an infinite number of operators as is necessary to allow 
a smooth joining to the WRGE solution. So doing, the full effective potential 
is then constructed as 

Veff{x, y) = V{x, y, f) + Vnoopix, y, f) , (10) 

where Vi^ioop{x,y,t*) is the analytically calculable 1-loop effective potential 
which accounts for the effects of the degrees of freedom carrying momentum 



8 



Figure 2: Ai (t) as a function of t computed using the nonperturbative WRGE 
with fixed Yukawa couphng. 

less than e^**A. A Coleman- Weinberg instability is signalled by a non-trivial 
global minimum of Veff{v) = Veff{x,y)\y^i^2^i^4 with vanishing renormal- 
ized mass. The condition of vanishing renormalized mass, V^^^\x=y=o-, is 
achieved by an appropriate tuning of the parameter m^(0) appearing in 
Eq. (2). If such a minimum appears v — e~*°A, then the system can 
sustain a hierarchy only over to e-foldings. 

Focusing on the specific initial couplings of Ai(0) = 0, A2(0) = 10, (?^(0) = 
1, we numerically integrated the WRGE and found that for t* ~ 1.5, V{x, y, t) 
was linear in t — t* and smoothly joined onto the 1-loop effective potential. 
As an indication of this, we plot in Figs. 2-4 the behavior of the coefficients, 

AiW = ^Vgf\^^y^oA2it) = ^Vy^^^.=y=o and = e-2*(A) V^/ZU-^.^o, 
of the two independent dimension four operators and one of the dimension 
six operators as a function of t. Note that we have included the explicit factor 
of ~ accompanying the dimension six operator in the definition of 
its couphng. As is clearly demonstrated, these exhibit a hnear behavior in 
t — t* for t* ~ 1.5. The behavior of these couplings in this region is identical 
to that obtained using the linearized in t — t* 1-loop effective potential with 
anomalous dimensions neglected. Note that for small t values, the coefficient 
of the induced dimension six operator (and other irrelevant operators) are 
quite sizeable and play an important role in the dynamics. As t increases, the 
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Figure 3: A2(t) as a function of t computed using the nonperturbative WRGE 
with fixed Yukawa couphng. 

effects of the irrelevant operators is diminished and the system is eventually 
attracted to the space spanned by the relevant and marginal operators only. 
Thus the contribution of the degrees of freedom with momentum less than 
e~**A can then be included using the 1-loop approximation. 

Since the WRGE approach includes an infinite number of operators, it 
would require an infinite dimensional space to plot the renormalization group 
flows. Thus we focus instead directly on the vacuum effective potential, 
Veff{v), defined in Eq. (10). In Fig. 5, we plot the generalized Yamagishi 
function Y(t) — — ^e"'*— ^ as a function oi t — —in v. A first order 
transition is signaled if this function crosses the stability line Y = 0. It is 
seen that such a zero occurs at ~ 2.0 and corresponds to a non-trivial global 
minimum of the effective potential at f ~ 0.14A. The value of to was seen 
to be insensitive to changes in the choice of t*. Thus the transition goes first 
order and a hierarchy of only ~ 2.0 e-foldings can be estabhshed. This result 
is in qualitative agreement with that found using the 1-loop perturbative 
approximation. 

We have also investigated the phase structure of the model for other values 
of the initial coupling parameter space using the WRGE with fixed g^{0). For 
instance, for Ai(0) = — 2, A2(0) = 10, g'^(O) = 1, the Yamagishi function Y{t) 
crosses the stability line at to — 1.9 indicating a first order transition after 
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Figure 4: i{t) as a function of t computed using the nonperturbative WRGE 
with fixed Yukawa couphng. 



Figure 5: Y{t) as a function of t for the initial parameters Ai(0) = 0, A2(0) = 
10, g^{0) = 1 computed using the nonperturbative WRGE with fixed Yukawa 
couphng. 
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1.9 e- foldings. Once again, this is in qualitative agreement with the 1-loop 
perturbative result of to = 1-0. Similarly, for Ai(0) = 2, A2(0) = 10, ^'^(0) = 1, 
we find the chiral symmetry breaking transition turns first order at to = 2.1 
while the 1-loop perturbative approximation gives to = 1-7. 

Using the nonperturbative WRGE, we have studied the stability of fine 
tuned hierarchies in a f/(2)^ x U{2)r chiral model for a range of initial cou- 
plings satisfying A2(0) >> 5'^(0) ~ |Ai(0)| ~ 1. In general, we found that 
the chiral symmetry phase transition turns first order in close proximity to 
the compositeness scale. This is in qualitative agreement with the results 
obtained using the 1-loop improved effective potential. In performing our 
analysis, we were forced to employ various truncations in order to render 
the numerical integrations tractable. Clearly the sensitivity of our results 
to these various approximations needs further scrutiny. Moreover, some of 
these truncations must be eliminated before we are able to probe the region 
of initial parameter space with larger Yukawa couplings. We hope to address 
these issues in future studies. 

We thank Sergei Klebnikov for many enjoyable useful discussions and 
Bijan Haeri for assistance with the numerical computations. This work was 
supported in part by the U.S. Department of Energy under grant DE-AC02- 
76ER01428 (Task B). 
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